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Abstract — We consider proportional fair rate allocation in a 
heterogeneous network with a mix of LTE and 802.11 cells which 
supports multipath and multihomed operation (simultaneous 
connection of a user device to multiple LTE BSs and 802.11 
APs). We show that the utility fair optimisation problem is non- 
convex but that a global optimum can be found by solving a 
sequence of convex optimisations in a distributed fashion. The 
result is a principled approach to offload from LTE to 802.11 
and for exploiting LTE/802.11 path diversity to meet user traffic 
demands. 


I. Introduction 

In this paper we consider rate allocation in a heterogeneous 
network with a mix of LTE and 802.11 cells. Integrated 
design of LTE and 802.11 is topical in view of the continuing 
increases in data traffic, the fact that many cellular network 
operators also operate a large network of 802.11 hotspots 
and that user handsets are now typically equipped with both 
LTE and 802.11 interfaces. Rather than an either/or handover¬ 
like setting where the question is which network to use our 
interest is instead in settings where user devices jointly use 
the LTE and 802.11 networks and may send data across both 
simultaneously. Eurther, we consider situations where user 
devices may connect to multiple LTE Base Stations (BSs) and 

802.11 Access Points (APs) simultaneously. This allows us to 
encompass the Coordinated Multi Point (CoMP) features in 
release 11 of LTE 0 which allow coordinated transmission 
and reception across multiple BSs and also multihoming to 
multiple APs and through a single WLAN card as proposed 
in, for example, Q. Of course connection to a single BS or AP 
remains as a special case within this more general framework. 
Simultaneous transmission across multiple interfaces might be 
implemented by routing each connection across one or other 
network (in a manner akin to load balancing) or by striping 
connections across both networks (e.g. via use of transport 
layer protocols such as multi-path TCP 0). 

Our focus in this paper is on how to allocate the avail¬ 
able LTE and 802.11 bandwidth amongst user devices in a 
heterogeneous network, and in particular how to determine a 
proportional fair rate allocation. Our main contributions are 
as follows. Eirstly, we develop a throughput model for het¬ 
erogeneous networks that include both LTE and 802.11 links 
and which encompasses multipath and multihomed operation 
(simultaneous connection of a user device to multiple LTE 
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BSs and 802.11 APs). We show that the rate region is non- 
convex, and is also not log-convex. Our second contribution 
is therefore a sequential convex optimisation approach, based 
on determining a sequence of maximal convex subsets, that 
we show converges to a global optimum. This optimisation 
approach is suited to distributed implementation. Thirdly, we 
present a number of application examples demonstrating how 
this framework can be used to develop principled approaches 
to offload from LTE to 802.11 and for exploiting LTE/802.11 
path diversity to meet user traffic demands. 


A. Motivating Example 

Consider the simple example in Eig [T] where the network 
consists of one cellular BS and one 802.11 AP. UE U 2 is 
located close to the 802.11 AP and so uses a physical rate 
of 54Mbps. UE ui is further from the AP and so uses a 
lower physical rate of 1Mbps. Both UEs are located a similar 
distance from the LTE base station and use the same physical 
rate of 10Mbps for their LTE connection. The physical rates 
available on the 802.11 and LTE links are summarised in Table 
[I] Eor simplicity, we assume both UEs are saturated i.e. always 
have a packet to send. 
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Eig. 1; Example of a mixed LTE/WiEi network, with one 

802.11 AP, one LTE BS and two UEs ui and U 2 each equipped 
with both an 802.11 and an LTE interface. 


When the UEs are contending for the available resources, 
the achievable throughput is determined by the channel ac¬ 
cess method. For 802.11 the random access approach leads 
to collisions and reduced throughput when more than one 
station contends for access. When ui and U 2 use only the 

802.11 WLAN, their data rates are, respectively 0.46Mbps and 
24.9Mbps for the standard 802.1 In MAC settings (recall the 
UE physical rates are 1Mbps and 54Mbps). In the case of 
LTE, time-frequency slots are allocated to UEs by the cell 
BS. When the UEs use only the LTE cell, the BS allocates a 
rate of 5 Mbps to each UE. 

It can be seen that ui enjoys better throughput via the LTE 
link than over 802.11, while U 2 ’s throughput over the 802.11 
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link is significantly better than that over the LTE connection. 
Due to this connection diversity the potential exists for more 
efficient decision-making to improve performance. For this 
example we later show that if ui sends all traffic via the 
BS, and U 2 sends all traffic via the AP, rates of 10 Mbps 
and 45 Mbps can be achieved by the UEs and that this is the 
proportional fair rate allocation. This is summarised in Table |I] 
Note that simply splitting the traffic for each UE equally over 
the 802.11 and LTE networks would yield rates of 5.46Mbps 
and 29.9Mbps i.e. worse for both users than the optimised 
allocation. 

TABLE I: Motivating Example: Achievable data rates 



PHY Rates 
[Mbps] 

Technology Rates 
[Mbps] 


BS 

AP 

LTE Only 

802.11 Only 

Optimised Multi-RAT 

Ui 

10 

1 

5 

0.46 

10 

U2 

10 

54 

5 

24.9 

42.4 


II. Related Work 

The problem of resource allocation in heterogeneous wire¬ 
less networks (HetNets) has been the subject of signihcant 
interest in the literature in the recent years. Much of this 
work has focussed on network selection, namely the task of 
selecting the best network for a UE to use. See Q for a 
survey on network selection in HetNets. In network selection 
two modes of UE operation are typically considered: multi¬ 
modal and multi-homed. Although both multi-mode and multi¬ 
home operations support multiple radio interfaces, only the 
latter supports multiple TCP flows across disjoint paths. The 
task of network selection is usually based on a specihed 
utility function with various methods (combinatorial, fuzzy 
logic, MDP, game theory, etc.) proposed to solve the resulting 
utility-based network selection problem. For example, 0,0 
consider combinatorial optimisation, 0 considers a linear 
utility function in a multi-mode UMTS/GSM network and 
0 a linear utility function which is the sum of logarithmic 
utilities of individual BS-UE pairs. In ||7) a network-user asso¬ 
ciation problem in a WLAN/UMTS hybrid network for multi- 
mode operations is solved using Markov Decision Processes 
(MDPs), although the complexity scales poorly with network 
size @.In0 a non-cooperative game is formed among users 
of two different classes corresponding to 3G/LTE and 802.11. 

Our work differs from this previous work in that (i) we 
consider concurrent association of the users among OFDMA 
LTE BSs and 802.11 CSMA APs by considering both access 
network models, (ii) we formulate the problem as a propor¬ 
tional fair resource allocation problem, (iii) we show that 
although the problem is non-convex it can be solved by solving 
a sequence of convex optimisation methods and (iv) we use a 
low-complexity distributed optimisation method to solve each 
convex optimisation problem. 

III. Network model 
A. Network topology 

We consider a fairly general multi-RAT setup where UEs 
can potentially connect to multiple LTE base stations and 


multiple 802.11 APs simultaneously. Networks where UEs 
are constrained to connect to a single LTE BS and/or 802.11 
AP are then special cases within this general setup. We note 
that the Coordinated Multi-Point (CoMP) features in release 
11 of LTE already allow transmission and reception across 
multiple BSs g, in. For 802.11 devices, connection to 
multiple APs might be implemented using the AP-hopping 
approach described in Q, where an 802.11 client periodically 
hops between a set of APs, informing those APs with which 
it is not currently transmitting/receiving that it is in power- 
save mode so that these APs buffer packets directed to the 
client until it next connects. This software implementation is 
a client-side only change that does not require more than a 
single 802.11 radio and which uses standard hardware. 

Let A denote the set of 802.11 access points and U := 
{1,2,... ,U} the set of users. Associated with each access 
point a S is the set of users Ua dU which can associate 
to it (so capturing geographic and access control constraints). 
We also let Au := [a G A : u € Ua\ U {0} denote the set 
of access points with which user u can associate, where 0 
corresponds to the case where the user is not associated with 
any access point. This also dehnes a bipartite graph Q with 
one set of vertices corresponding to users, a second set of 
vertices corresponding to the access points and edges between 
each client vertex and the set of access point vertices to which 
it can connect. 

Similarly, let B denote the set of LTE base stations, Uh the 
set of users located in the cell operated hy b G B and Bu the 
set of base stations with which user u can connect. 


B. LTE throughput 

Let I denote the set of available LTE sub-channels and 
let Q ^ denote fraction of time sub-channel i of BS 6 is 
used by user u. Let be the set of user-BS pairs for which 
transmissions interfere on sub-channel i (this dehnes a conhict 
graph). We consider joint sub channel allocation and time 
sharing so that at a given time interfering BSs do not transmit 
on the same sub-channel GD- Hence, the total allocation for 
each sub-channel must satisfy the following: 

5] Cm<i (1) 

(u,b)GE‘^ 

When there is no frequency reuse, this constraint simplihes 
to Eueu u — Vf S X. Letting w* denote the sub¬ 

channel bandwidth, then the achievable rate of user u is given 
by: 


beB iex 


UJ log(l-b —) 


yuGU ( 2 ) 


with ylu = denoting the SNR on sub-channel i of 

BS b where p\, e\ „, cr„ denote the BS power on channel 
i, channel gain, and noise power at the receiver respectively. 
Parameters jSi and /?2 reflect the LTE bandwidth and SINR 


implementation efficiencies 112|. Hence the rate region of the 
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LTE network is the following 

Tlite = |r = XI XI Cb,ul^i^" log (1 + 

bGB iGl 

r<ru<f, 0 < Cb,u < 1: 

XX^M<l>VrGl| (3) 

u^Ub^B -1 

where r G is the vector formed by stacking the user 
throughputs ru, u GU and we have also added the constraint 
that each user has maximum and minimum rates f and r 
respectively. Note that r may be 0 and f may be oo. It can 
be seen that the LTE rate region is convex (a polytope in fact, 
since the constraints are linear) in the time-frequency sharing 
factors du- 

We note that this model can be extended to include dynamic 
power allocation per sub-channel, for example using a similar 
approach to E) and will still be convex. 


C. 802.11 WLAN scheduling 

We consider an AP hopping approach where time is par¬ 
titioned into scheduling slots of duration T indexed by f = 
1,2,... such that in slot t user u G U operates its 802.11 
MAC in association with at most one access point (reflecting 
the constraint that user u only has a single 802.11 radio). 
Extension to allow simultaneous connection to multiple access 
points is straightforward. 

The access point which user u selects at time slot f is a 
random variable A„ t which takes value a from Au \ {0} 
with probability Za,u (so 0 < Za,u < 1 and < !)■ 

Access points are selected independently at each slot and by 
each user. Note that with this randomised schedule a client 
is not associated with any access point in a time slot (i.e. 
A„,t = 0) with probability := 1 - 


set S which have cardinality k. Therefore Vn-i{lda \ {u}) is 
all subsets of Ua\{u} with cardinality n — 1. Define 

Pa,u,n ■= ProbCNa.t = n\Au,t = a) (4) 

to be the probability of n users being associated with AP a 
conditioned on user u being associated with a. We have 


n (1 - Za,v) n = 1 

Pa,u,n = ■( n Za,v H (1 “ ^a,v) n > 1 

Ua€ vGUa V^Ua\{l^a ,'i^} 

\.Vr,-l(Ua\{u}) 

(5) 

Letting ■> rewritten equivalently as 

1 


Pa^u.l — 


n ii+Wa,v) 


(6) 


and for n > 1 


Pa,u,n — 


X n 


Y[ (1 - Za,v) (7) 


_ 1 — Zn V 

«ae VeUa ’ ^eUa\{u} 

E n 'Wa,v 

Wa6'P„-l(Wa\{«}) VeUc 

n {^+Wa,v) 


n = 1 

E n '^a,v n> I 

UaGVr^-l{Ua\{u}) vGUa. 


(8) 



(9) 


then the expression for probability Pa,u,n simplifies to 

Pa.,u,n — Qa,u,nPa,u,lj ^ — 1, 2, . . . (10) 

Note that El!il Pa.«.n = Pa.«,l El!il = 1 “d SO 

ylWal ^ 1 

Z^n=l 'ia,u,n ^ ^ ■ 


t=l t = 2 t = 3 

Al,t 

A2,t 

Ni,t 

N2,t 

Pig. 2: Illustrating the 802.11 random variables A„ j and 
Na t related to user-AP association, for a two user, two AP 
example. 

Let random variable Na,t denote the number of clients 
associated with access point a in time slot t. This is illustrated 
in Fig for a setup with two UEs, and two APs oi, 02 . In 
slot 1 both UEs are associated with AP ai. In slot 2 user ui 
is associated with AP 02 , while user U 2 is associated with AP 
oi and in slot 3 both users are associated with AP 02 . 

By construction, random variables Na,t, f = 1,2,... are 
i.i.d. Let 2*^ denote the set consisting of the subsets of set S, 
and let Vk.{S) = {s : s G 2*^, |s| = k} denote the subsets of 


ai 

a2 

a2 


ai 

fli 

02 


2 

1 

0 


0 

1 
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D. 802.11 MAC slots 

The number of users actively using an AP varies from time 
slot to time slot. We therefore have to take some care to define 
the MAC slots within each WLAN appropriately. Within each 
WLAN MAC slots are induced by carrier-sensing, which 
may be idle slots of duration cr (no user transmits) or busy 
slots of duration Tf, (at least one user transmits, this includes 
both successful transmission and collision times). In order 
to distinguish between MAC slots associated with different 
access points, we index MAC slots by i G A x N i-e. the MAC 
slots associated with access point a are i = (a, 1), (a, 2),.... 

We let A4a,t denote the set of MAC slots of access point a 
that are fully contained within time slot t and be the set 
of MAC slots which are only partially contained in time slot t. 
Since the MAC slot duration is random, generally there will 
be boundary MAC slots which are only partially contained 
in time slot t but there are at most two such slots. We will 
assume that the time slot duration T is sufficiently large that 
partial MAC slots A4a,t can be neglected when calculating 
user airtime and throughput. This is illustrated in Figure 
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Fig. 3: Illustrating 802.11 time slotting : each user u can hop 
between multiple APs, but is associated to the same AP for 
time slots of duration T. Within a time slot of duration T, 
802.11 carrier sense defines MAC slots. Ala.t denotes the set 
of MAC slots fully contained within time slot t. 


E. 802.11 throughput 

Let La,u, Tsucc,a,u, TcoU and a denote, respectively, the 
payload in bits of a successful transmission by station u 
when associated with AP o, the mean duration of a successful 
transmission, the mean duration of a collision and the idle 
slot duration. To simplify notation we confine consideration to 
client uplink transmissions. Extension to include AP downlink 
transmissions is straightforward, simply requiring appropriate 
book-keeping of which fraction of AP transmissions is directed 
to each client. For simplicity we also assume that Tsucc,a,u 
is the same for all users transmitting in WLAN a and let 
T — T 

^ succ,a — ^ succ.,a,U‘ 

In time slot t let Ua,t C Ua denote the set of users associated 
with access point a and recall that = \0{a,t\- Let „ be 
a random variable which is equal to 1 if station u transmits 
in MAC slot i G Aia,t and 0 otherwise where u G Ua- We let 
Xi,u — 0 for users u ^ Ua^t not associated with access point 
a in slot i G Ma.t- For u G Ua,t we assume that the {Xi^u} 
are i.i.d. in both u and i i.e. Xi^u ^ X, and we let r := 
Prob{X = 1). This can be achieved in practice by setting 
CWmax = CWmin and using the same value of CWmin for 
all associated users, in which case r = 2/(C'Wmm + l)- We let 
yi^u be a random variable that equals 1 if station u transmits 
successfully in MAC slot i, and 0 otherwise. We let Zi equal 
0 when one or more users transmits in MAC slot i, and 1 
otherwise. Since transmissions {Xi^u} are i.i.d. for u G Ua,u 
i e we have that ~ with = 

t(1 - and Z, ~ with = (1 - 

The throughput of user u in WLAN a is given by 


Sa,u = lim 

k—¥tx> 

= lim — 
k—^oo k 




kT 


E E 


Ln 


tG{sG{l,2,...,fc}:A„,a=a} iGMa.t 

\Ua\ ,j-k I , j 

= lim y ^ 1 y y :y;,_ 

k—^oo k 

n=l (gT'' ieMa t 


T 


( 11 ) 

( 12 ) 

(13) 


Defining 


n 


Psucc,n = E^[^"] = 

(14) 

P,die,u = EiZ'-] = (1 - r)” 

(15) 

Pcoll,n — 1 Pidle,n Psucc,n 

(16) 

have 


lim 'y yi u 

(17) 

= E[\Ma,t 1 |N,.t = n, = a]E[y’-] ^ 

(18) 

Pidle.,n^ P Psucc,nP'succ,a P Pcoll,nPcoll 

(19) 


\j-k I 

We also have limfc_j.oo = Prob{Au,t = a, Na,t = n) = 
Prob{Au,t = a)Prob(Na,t = ri\Au,t = a) = Za,uPa,u.n- 

Thus, the throughput of user u in the WLAN operated by AP 
a is given by 


IWal 

'Lh 

n = l 


Pa,u,nT{l - tY" 


idle,n 


O' + Psi 


+ Pc 


coll,nX coll 


( 20 ) 


1) Each User Associates to a Single WLAN: It is helpful 
to consider the special case where each user u connects only 
to a single WLAN a, in which case „ = 1 and Za\u = 0 
for a' Y Letting no denote the number of users connected 
to WLAN a then Pa,u,na = 1 and Pa,u,n = 0 for n Y kia- L 


follows that (201 simplifies to. 


-(1-t)”“-1Lo.. 


PsucCjTia^succ^a Pcoll^nJ^coll 


( 21 ) 


which is identical to classical 802.11 WLAN models e.g. that 


in 113|. However, when the number of users using a WLAN 
varies from time slot to time slot, (|^ yields the shared 
throughput. 

2) Change of variables: It will prove useful to work in 
terms of the quantities f = y: and Wa^u = rather 

than r and Za.u^ respectively. Also we will use the normalised 

■ g 

throughput pa.u = where Ca^u = Recalling 

Pa,u,n = <la,u,nPa,u,i and letting Na = we have 


Pa,v 


1 + Wa.u 


\Ua\ 

E- 


Qa,u,r 


f Ts, 


n (1 + Wa,«) <!>«(<)>) Tc 


( 22 ) 


with 


where := {s S {1, 2,..., fc} : Au,s = a, N,., = n} and ^ ^ ^ + i)« _ i (23) 

we have used the fact that „ = 0 when u ^Ua- 
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3) 802.11 Rate Region: With the above change of variables 
the rate region of the 802.11 network is, 

T^wifi — I ® ^ ^ Pa^u^a.u-) 

IWal 


V£Ua 

Qa,u,n — E n Wa,v, n> 1, 


Wae veUc, 

•p„_i(WA{«}) 


E 

aeA 


1 + Wa 


^ I 5 0 ^ Pa,u — I 5 ‘^a.u ^ 0; Qa.u^n ^ 0 


Vn gU, Wa G a, \/u gU 


(24) 



(a) 1 AP: Tai = 0.95 (b) 2 APs: Ta^ = 0.95, Ta^ = 0.2 


Fig. 4: Examples of 802.11 rate region IZ^^ifi where hopping 
between APs is admissible. In (a) two users hop between 
connecting and not connecting to an AP and in (b) two users 
hop between two APS for which the transmission probabilities 
differ. 


where s G is the vector obtained by stacking the user 
throughputs s^, u G U. Letting pa,« = logpa.^, Wa.u = 
logWa.u, we also have that the rate region of the 802.11 
network can be written as 



■ — 


,U5 

1 


a^Au 





PA 

Pa,u E 

-E 

log(l+e”-'’ 

0 + log^ 


V&Aa 


n—1 

Qa,u,n — 

E 

n 

n > 1 


ZkA 

T, $„ 


Mae UGWa 


E 

q-GA 




1+6-“ 


— Qa,u,n — 0; 


Vn Ghl, Va e A, GU 


(25) 


4) Example 802.11 Rate Region : Recall the simple exam¬ 
ple in Fig where we have one BS, one AP and two users. 
The 802.11 rate region is, 


In the special case (see above) where for each user u we 
force Za,u = 1 for one AP a, then each user associates to 
a single, fixed AP and the only design freedom lies in the 
LTE allocation which is a convex problem. However, when we 
additionally allow the UEs to adjust the fraction of time that 
they connect to this single AP then the rate region becomes 
non-convex, as illustrated in Figure Qa), and when users may 
connect to more than one AP is illustrated in Fig Qb). Fig 
Qb) plots the rate region for networks with two APs 

and two UEs. The non-convexity of the rate region in Figj^b) 
is evident. Note that the two extreme points of the “notch” in 
the rate region indicated in the figure correspond to parameter 
values = [1 0 0 1 ]^, and z^ = [0 1 1 0 ]^, where z := 
[•^ 1,1 -^ 1,2 -^ 2,1 •^ 2 , 2 ]^- That is, point 2 } corresponds to ui being 
connected to AP oi 100% of the time and U 2 to AP 02 , while 
z^ is the reverse. The UE rates along the time-sharing chord 
aA -f (1 — a)z^, a G [0,1] are indicated by the line marked 
on Fig Qb). It can be seen that this lies in the interior of 
the rate region. Simple time-sharing therefore does not yield 
convexity, and this is due to the collisions in 802.11 WLANs 
shared by more than one UE. 


R-Wifi — j S : Si — S2 — Pl,2Cl,2 

^ Wl,l / 4> , 

- (l-bwi,i)(l-f wi.2)^4-i T, ’ 

91 , 1,2 < Ull, 2 , 


Pi,2 < 


Wl,2 ( ^ \ ^ 'iTsi,2 

( 14 - -b wiy)' ^ 


91,2,2 < 


0 < Pi,u < 1, Uii.u > 0, 91 , 1 *,n > 0, {u,n|u,n G {1,2}}| 


( 26 ) 


F. Non-convexity of Network Rate Region 

The network rate region is the joint LTE/802.11 rate region 
IZite X TZwifi- The LTE rate region IZue is convex and so the 
convexity/non-convexity of the network rate region depends 
on the convexity/non-convexity of the 802.11 component of 
the network rate region IZwifi- 


IV. Proportional Fair Rate Allocation 
A. Utility Fair Optimisation 

We are now in a position to consider the main focus of this 
work, namely finding a proportional fair rate allocation for the 
joint LTE/802.11 network. This is the solution to the following 
utility optimisation P, 

max V log {su+Tu) (27) 

s,r ^ 
uGU 

s.t. s GUyjifi, r GlZite (28) 

As already noted, the rate region IZwifi (and so IZwifi) 
is generally non-convex and so the utility fair optimisation 
P is non-convex. However, as we will see the optimisation 
possesses sufficient convex structure to allow near-optimal 
solutions to be found in an efficient manner. We proceed by 
first considering methods for approximating a non-convex set 
by a maximal convex subset. In this way we can define a 
convex approximation to P for which solutions can be effi¬ 
ciently found. We then consider a concave-convex procedure 
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Fig. 5: Illustrating a maximal convex subset (indicated in 
darker shade) of concave set y < e^. 


that adaptively selects the maximal convex subsets with the 
aim of maximising the optimisation objective. This procedure 
is guaranteed to converge to a stationary point of non-convex 
optimisation P. Although this stationary will in general be 
sub-optimal, in practice we hnd that it is usually near-optimal. 


B. Approximate Optimisation Via Maximal Convex Subsets 

1) Maximal convex subsets: Sets of the form i? := {(x,y) : 
y < G^,{x,y) G K^} will prove important in our analysis. 
Set E is concave, see Figure Nevertheless, large convex 
subsets of this set can be readily identified. In particular, the 
set F := {(a:, y) '■ y < e^{l -\- x — x), {y, x) € K^}, with 
parameter a; S M, is convex with F C E, as illustrated in 
Figure Convexity follows from the fact that the constraint 
dehning set F is linear in x and y. That F is contained in 
E follows from the observation that the complement E := 
{{x,y) : y > e^,{x,y) G of set F is convex and the 
boundary of F" is a supporting hyperplane to F at point (x, e^). 

2) Maximal Convex Subsets of 802.11 Rate Region: The 
main constraints in the 802.11 rate region itmfi can be written 
in standard form as 

6 \ 

Pa,u - Wa,u + X! X! < 0 

v&Aa / 


Qa,u,n 


E n 

MoE VeUa 

■P„-l(Wa\{«}) 


< 0 


(29) 

(30) 


The terms on the LHS of the pa^u constraint are all convex. 
However, the second term on the LHS of the qa,u,n constraint 
is concave. We proceed by hnding a maximal convex subset. 
Letting 

^ • {(^a.u.nj \.^a,v } vCiVia ) ' ^a,v € M, C Ha j 

qa,u,n ^ E n e’"”’"} (31) 

iiaE t)GWa 

■p„_i(WA{«}) 

then a maximal convex subset is 


F' ■ — { (^a,u,n 1 }vGlta ) • ^a,v G 1^,17 G Ha ^ 

qa,u,n < ^ ^Ua,v 

WaE V&Aa. 

■p„_i(W„\{«}) 


(32) 


where „ := {lFY.v^Ua.\{u}i.'^cL,v-Wa,v)) and Wa,u 

is a design parameter. 

The 802.11 rate region also contains the constraint 

- E < 0 (33) 

It can be seen that the LHS is concave. Once again, we adopt 
a maximal convex subset approach 

^ ^ ' (1 -f pa.u Pa^u)t:a^u — 0 

aeAu 


where „ is a design parameter. The LHS is now convex in 
Pa,u as required. 

Lastly, the 802.11 rate region contains the constraint 


E 

aeA 




1 -be’' 


< 1 


(34) 


the LHS of which is sigmoidal in form, so neither convex not 
concave. Rewriting equivalently as 

observe that the first term in the RHS is convex while the 
second term is concave. Linearising the second term yields 
the following inequality 



a£A 


/ 2e2®“.“ 


p3Wa.,-. 


(1 + 


{Wa,u - Wa,u) < 1 

(35) 


which defines a maximal convex subset. 

3) Solving the Convex Optimisation: Using the maximal 
convex subset approach described above we obtain a convex 
optimisation the solution of which is feasible for non-convex 
optimisation P but will, in general, be sub-optimal. The 
Slater condition is satisfied and so strong duality holds. We 
can therefore solve this convex optimisation in a distributed 
manner using a primal-dual subgradient approach. 

Let X := G K" be the vector 

obtained by stacking LTE and 802.11 rate region variables with 
^ g ]a|i|x|B|xC/^ ^ ^ g ]^\A\xU^ q g ]^|^|xc/x(7 denoting, 

respectively, the vectors with elements Q Pa,u, Wa,u, qa,u,n 
fox h G B, i G X, a G A, n G H, u G H. We re-write the 
optimisation problem P in the following form, 

min /(x) (36) 

X 

s.t. /i^’^(x) — ( 7 ^*^(x) < 0, * = 1,2, ...,Z (37) 

with / : M" —>• K, : M" —>• K, : K —>• K being convex 
continuously differentiable functions, expressions for which 
are given in Appendix A. Let —p(9(x;x) be the maximal 
convex subset expression for non-convex function — 5 *^*^(x). 

(x; x) = (x) - (x) (x - x) (38) 

Then the approximate optimisation problem P^ is given by 
min /(x) (39) 

X 

s.t. /i^’^(x) — p^’^(x;x) < 0, *=1,2,...,/ (40) 
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Letting A := • • • , denote the set of multipliers as¬ 

sociated with the rate region constraints 1 to I, the Lagrangian 
for optimisation problem is 

i 

L(x,A;x) = /(x)+^A«(/rW(x)-pW(x;x)) 

i=l 

The standard primal-dual subgradient approach in Algorithm 
[1] can then be used, for example, to hnd a solution to 
optimisation P^- 


Algorithm 1 Distributed primal-dual algortithm 

Initialise: t = 0, x(0), A(0), step size a > 0 

do 

x(t + 1) = x(t) — adxL(x (t) , A (t) ; x) 

A(t + 1) = [A(t) + adAL{x (t) , A (t) ;x)]^ 
t i — t -|- 1 

loop 


4) Message Passing Required: 

1) To update the LTE sub-channel airtime C,l ^ each UE 
needs (i) the SNR to its LTE BSs (which it already 
knows) and (ii) the multiplier associated with the sub¬ 
band constraint (which can be communicated by a BS). 

2) The 802.11 association probability Wa,u and WLAN 
parameters qa,u,n can be updated using information avail¬ 
able locally at AP a (no need for message passing) 
together with knowledge of the multiplier associated with 
the constraint that association probabilities for user u 
sum to one. This requires that Wa,u as well as Wa,u be 
communicated by AP a to UE u and the multiplier then 
communicated back from UE u to AP a. 

3) The 802.11 normalised throughput „ can be updated 
by UE u using local information together with knowledge 
of the multiplier for the rate constraint on (which can 
be communicated by AP a). 

C. Adaptation of Maximal Convex Subsets 

The convex optimisation in the preceding section yields 
a solution which is feasible for non-convex optimisation P 
but which is, in general, sub-optimal for P. The degree of 
sub-optimality is dependent on the choice of maximal convex 
subsets used to derive the convex optimisation, e.g. when the 
convex subsets contain at least one point which is an optimum 
of the non-convex problem P then the solution to the convex 
optimisation P^ will in fact be optimal for P. 

When we have some knowledge of the part of the network 
rate region where the optima of problem P are likely to 
lie, we can use this to select the maximal convex subsets. 
Such knowledge might be available from prior experience, e.g. 
previous solutions to similar network conhgurations, or from 
structural insight. In the next section we illustrate how such 
prior information can be used to obtain near-optimal solutions. 

However, it is also possible to automate this process and 
this is the focus of the present section. The basic idea is to 
iteratively update the choice of maximal convex subsets based 
on the current solution to the convex optimisation. Having 
obtained the solution to the convex optimisation for the current 


choice of maximal convex subsets, one natural approach is to 
use the components Wa,u, Pa,u of this solution as the values 
for the parameters „ and in this way to dehne a new 

maximal convex subset. 

1) Convergence: In more detail, let C := {x G B : 
— pb)(x) < 0, i = I,-- - .,^} with B C K" non¬ 
empty, convex and compact (closed and bounded). As we will 
see, set B is needed for technical reasons, to ensure that C is 
compact, but B can otherwise be chosen arbitrarily and can 
be viewed as augmenting the set of convex constraints 
Of course, we assume that set C is non-empty. Now consider 
the iterative update 

Xfe+r G D^^f{x), k = !,■■■ (41) 


with xi G C, 




argminxec., /(x) if /(x^) = f*^ 

{x G Cxfc : /(x) < /(xfe)} otherwise 

(42) 


and f*^ := minxec^^^ /(x), Cx^ := {x G B : /ib)(x) - 

5«(x;xfe) <0}. 

Observe that the solution to convex optimisation Pxt lies in 


set iAxfc/(x) and so can be used in update (41 1 . However, set 


iAxfc/(x) dehnes a wider class of updates, including suitable 
approximate solutions to P^^ that may be easier/faster to hnd, 
and generalises the concave-convex approach of M- 
We have the following convergence result. 

Theorem 1 (Local Convergence): Iterative update ( |4T] l con¬ 
verges to a stationary point of non-convex optimisation P. 

Proof: See Appendix B. ■ 


By adding a stochastic search component to update (41 1 we 


can strengthen this to obtain a global convergence result. 

Theorem 2 (Global Optimisation): Consider iterative up¬ 
date 


Xfc+i G 


‘-fe+r 


w.p. e 


Pxfc / (x) Otherwise 


Vk+i = argmin{/(?/fe),/(xfe+i)} 


(43) 


with Xi G C and where {X^} is a sequence of random 
variables taking values in C satisfying Pro5(Xfc_|_i G E|Xfc G 
T) > pp(T n E) for any T, E C C where p{T n E) denotes 
the volume of their intersection and parameter rj > 0. Then, yk 
converges to the optimum of optimisation P with probability 
one. 

Proof: See Appendix B. ■ 

Update ( [43] l is no longer purely a descent update, but rather 
with probability e an update is made which may lead to the 
objective / increasing, allowing escape from unfavourable 
stationary points (saddle points etc) and from local minima. 
The requirement on random variable X^+i can be satished 
by, for example, selecting X^+i uniformly at random within 
a ball about the current point Xj,. 

Although the proof of Theorem could be used to upper 
bound the convergence, this bound would be very loose and 
so while Theorem provides some comfort as to the ability 
of update (43 1 to hnd a global optimum, evaluation of the 
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Fig. 6: Example illustrating global convergence, a = 0.1, e = 

0 . 2 . 


convergence rate really needs to be carried out via numerical 
experiments. 

Note that keeping track of the best Xk+i at each step of 
Update ( |43| l doesn’t lend itself to distributed implementation 
and so we drop that step in the proceeding examples (although 
this also means that we lose the convergence result of theorem 
2 ). 

2) Descent Step: Descent step can be realised 

in many ways. One is by use of a truncated version of the 
primal-dual subgradient approach: 


Algorithm 2 £>xfc/(x): Truncated Primal-Dual Updates 


do 


z(i) = x(t) 

z{t + 1) = z(t + 1) — ad-^L(z (t) , A (t) ; xj.) 
aO) {t + l)= [aO) (i) + adAL{z (t) , A (t); x^)] ^ 
t = i + 1 


until f(z{t + 1)) < /(xfc) and z{t -I- 1) S Cxj, 
XJ. + 1 = z(t + 1) 


While this primal-dual approach might not lead to a feasible 
point at every step, the requirement that x^+i be feasible can 
be relaxed to the requirement that (i) x^+i is bounded and 
(ii) Xfc+i is feasible for sufficiently large k. Theorem [T] then 
continues to hold, with only minor changes to the proof 
3) Example: The following simple example illustrates the 
convergence of update (j4^ to a global optimum. Consider the 
optimisation problem 

2 4 

min X — X 
xG[-4,4] 

It can be verified that in the interval [—4,4] this has a global 
minimum at x* = -EXjsfi with /(a:*) = —0.25, but it also 
has a stationary point at a: = 0. Fig [^a) compares update ( |43| l 
with update (|4^ when starting from initial condition a; = 0 
(a stationary point). It can be seen that update ( |4l] l gets stuck 
at this stationary point whereas update ( [4^ is able to escape 
and find a global minimum. Fig |^b) shows realisations of 
update ( |43] l for a range of initial conditions, illustrating its 
insensitivity to the choice of initial condition. 

V. Motivating Example Revisited 
Armed with our proportional fair problem formulation and 
non-convex optimisation tools we now revisit the example in 


Fig in more detail. 

A. Network Setup 

The network consists of one LTE base-station, one 802.11 
AP and two UEs. We assume that the AP uses 802.1 In fTS) 
with the settings detailed in Table and the A-MPDU frame 
structure illustrated in Fig For simplicity we assume that 
only two LTE sub-channels are available, denoted sub-band 
one and sub-band two. We consider LTE Frequency Division 
Duplex (FDD) systems with lOMHz system bandwidth and 
uplink transmissions with one data stream from single antenna 
UEs. The LTE system parameters are also detailed in Table 

im 


PSDU 

•s- 

PHY Header | A-MPDU 



I Subframe 1 | Subframe 2 | 



Delimiter 

MPDU 

Header 



MSDU 

FSC 

Padding 



I Subframe K | 


Fig. 7: 802.1 In frame format for Aggregated-Mac Packet Data 
Units (A-MPDUs). 


TABLE II: Network parameters. 


802.11 Parameters 

CWmin 

Physical Rate, {ca,i,Ca.2) 

Physical Rate, (ct,^i,Ch^2) 

Basic Rate (Rtasic) 

Slot Time (a) 

PLCP Preamble 

PLCP Header 

MAC Header 

Payload size (L) 

Number of MPDU subframes for u\ {Ki) 

Number of MPDU subframes for U 2 {K2) 

FSC 

Delimiter 

Padding 

ACK 

RTS 

CTS 

15 

(1, 54) Mbps 
(5,1) Mbps 

1 Mbps 

9p.s 

16fis 

48 bits ®Ca,u 

192 bits @Ca.u 

1500 bytes 

1 

54 

32 bits ®Ca,u 

32 bits @Ca,u 

32 bits @Ca,u 

112 bits @Rbasic 
160 bits @Rbasic 
112 bits @Rbasic 

LTE Parameters 

Duplex Mode 

System Bandwidth 

Carrier Frequency 

Number of Antennas at BS 

Number of Antennas at UE 

Guard Band Overhead 

DMRS Overhead 

Random Access Overhead 

Cyclic Prefix Overhead 

Bandwidth Efficiency (/3i) 

SNR Efficiency (/32) 

SNR ,(7‘*J,7S*^), i = l,2 

FDD 

2 X 10 MHz 

2 GHz 

1 

1 

10% 

14.3% 

0.625% 

6.66% 

71% 

1 

(4.9.4.9) dB 

Optimisation: 

Maximal Convex Subset 1 (xi) with -((Pa,!,pa. 2 ), 

{(-1.0).(-2,2)} 


Maxima] Convex Subset 2 (X 2 ) with | (pa,i) pa. 2 )i(fiia,i> 1 ^ 0 , 2 )} 

{(0,0) .(0,0)} 


Maximal Convex Subset 3 (X3) with {(pa,i,pa. 2 )>(ti'a.i.'*r'a, 2 )} 

{(0,0),(1,1)} 


Q (step size) 

0.01 


e 

0.2 


B. Rate Allocations 

Using the maximal convex subset Xi given in Table 
Algorithm [^yields rate allocations of 10 Mbps and 42.4Mbps 
for UEs ui and U 2 respectively. Using update ( |4T] l to adapt 
the maximal convex subset the rate allocation improves to be 
10 Mbps for UE ui and 46.5Mbps for UE U 2 - Details of the 
solutions found using the two approaches are given in Table 

m 

For comparison, when ui and U 2 use only the 802.11 
WLAN their data rates are, respectively 0.46Mbps and 
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24.9Mbps while when the UEs use only the LTE the BS 
allocates a rate of 5Mbps to each UE. Splitting the traffic for 
each UE equally over the 802.11 and LTE networks would 
yield rates of 5.46Mbps and 29.9Mbps. 


TABLE 111: Solutions given by Algorithm and update (41 1 . 


Pai'ameter 

Algorithm 1 Convex subset xi 

Update 441| 

Ui 

U2 

Ui 

U2 

Cb,u 

1 

0 

1 

0 

pa,u 

0.05 

0.85 

0 

0.98 

Tu 

10 

46.5 

10 

48.8 

Pa,u,l 

0.06 

0.98 

0.01 

0.99 

Pa,u,2 

0.85 

0.01 

0.97 

0.03 

^a,u 

0 

0.93 

0.01 

0.99 


C. Convergence 

Eigure illustrates the convergence of update ( [4T| l for 
various choices of initial convex subset, detailed in Table ||^ 
The first choice xi makes use of knowledge of the network 
to estimate the rate region in which the optimum is likely to 
lie. The second and third choices X 2 and X 3 are randomly 
selected. 

Eor comparison, Eigure shows the convergence of both 
updates ( |4T] l and ( |4T| l. Observe that in this example the extra 
complexity of update (43 1 does not yield a better optimum, and 
indeed we also observe this in the other examples presented 
below (we also carried out further tests, not shown here, that 
exhibit similar behaviour). This suggests that for the class of 
optimisation problems considered here update ( pT) tends to 
converge to a near optimal solution. 


VI. Example Scenarios 
A. 802.11 Offload 

Our problem formulation can be readily extended to provide 
a principled approach for offload of data trafiic from LTE to 
802.11. In an offload setting we would like to preferentially 
use the 802.11 network where possible. We can capture 
this requirement by augmenting the proportional fair utility 
function (^u + r„) with a cost related to LTE usage. 

Eor example, the cost function 


^(Cm) = EEE 


a 


b,u 


(44) 


GW bGB iGl /3lW* log (1 + -^) 


associates a cost with the LTE airtime used by each UE (the 
airtime being inversely proportional to the data rate). 

Consider a simple network setup with one BS, one AP and 
two UEs. Suppose the 802.11 physical rate is 54Mbps for both 
users and the LTE rates are 4Mpbs (corresponding to an SNR 


of 0.75dB), see Table IV Other 802.11 LTE parameters are 
detailed in Table |II] 

Using update ( [4l] i to solve the optimisation problem. Table 
[Vl summarises the solution found. The UEs share the 802.11 
AP and do not send traffic via the LTE network (the LTE rates 
fui = 0 = r^j). The throughput of each UE is 16.5Mbps, 
which is higher than the data rate of 2Mbps provided by the 


LTE network. Eigure 10(a) plots how the fraction Q b u 



(a) Throughput of UE ui 


(b) Throughput of UE U2 



Iteration 


I 

5 

= 

I 

I 

I 



(c) 2 a,l 


(d) 2a,2 


i 

I 

p 


0.8 

0.6 

0.4 

0.2 

0 


Convex Subset 1/Sub-band 1 —a— 
Convex Subset 2/ Sub-band 1 — e— 
Convex Subset 3/ Sub-band 1 —‘— 
Convex Subset 1/ Sub-band 2 — ■— 
Convex Subset 2/ Sub-band 2 • 

Convex Subset 3/ Sub-band 2 — ■— 


1 10 19 28 37 46 



(e)qiA = l,2 


(Oa 


1,2 


Eig. 8 : Convergence of update for various choices of 
initial convex subset, detailed in Table 



(a) Throughput of UE ui 




(b) Throughput of UE U2 



Eig. 9: Convergence of updates (41 1 and (431 when is 

drawn uniformly random within the ball of radius 0.025. 


LTE airtime used by each UE changes as the LTE data rate is 
increased, while Eigure [T0{b) shows the aggregate throughput 
(LTE plus 802.11) of each UE. It can be seen that for LTE data 
rates less than 5Mbps the LTE network is not used and data 
is fully offloaded to the 802.11 network. However, as the LTE 
data rate increases the LTE network is increasingly used to 
enhance the user throughputs. Observe also from Figure [T0|a) 
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TABLE IV: 802.11 offload example data rates. TABLE VI: 802.11 multihoming example data rates. 



PHY Rates 
[Mbps] 

Technology Rates 
[Mbps] 


BS 

API 

LTE only 

802.11 only 

Optimised Multi-RAT 

Ui 

4 

54 

2 

16.5 

16.5 

U2 

4 

54 

2 

16.5 

16.5 


TABLE V: Rate allocation for 802.11 offload example. 


Parameter 

Ui 

U2 

Tu 

0 

0 

Pa,u 

0.48 

0.48 

Su 

16.47 

16.47 

Pa,u,l 

0.046 

0.046 

Pa,u,2 

0.954 

0.954 

Za,u 

0.95 

0.95 


that the probability Za^u of associating to the 802.11 network 
remains almost constant as the use of the 802.11 network 
always results in an increase in the aggregate throughput. 



PHY Rates 
[Mbps] 

Technology Rates 
[Mbps] 


BS 

AP ai 

AP 02 

LTE 

WLAN ai 

WLAN a2 

Optimised Multi-RAT 

UI 

10 

36 

36 

5 

5.9 

8.9 

22.39 

U2 

10 

54 

18 

5 

8.87 

4.3 

22.39 


TABLE VII: Rate allocation for 802.11 multihoming example. 


Parameter 

Ui 

U2 

Cb,u 

0.5 

0.5 

po-l ,U 

0.0045 

0.92 

Pa2,u 

0.94 

0.0037 

Tu 

5 

5 

Su 

17.39 

17.4 

Pa\ ,u,l 

0.007 

1 

Pa2,u,l 

1 

0.007 

^ai ,u 

0 

0.99 

^0-2 ,u 

0.99 

0 


Once the 802.11 rate increases further, the extra capacity is 
then used by UE ui. 


1 

0.8 

0.6 

0.4 

0.2 

0 






(a) (b) 

Eig. 10: Illustrating 802.11 offload as LTE data rate is varied, 
(a): LTE airtime, (b): Aggregate LTE plus 802.11 UE through¬ 
put. 


(a) UE Ml (b) UE U2 

Eig. 11: Illustrating 802.11 offload as 802.11 rate is varied. 


We can readily extend consideration to include situations 
where UEs have a specified traffic load by adding additional 
constraints to the optimisation problem. Eor example, suppose 
that both UEs have a traffic load of 10Mbps video, the LTE 
available rate for each UE is 10Mbps (i.e. the total LTE 
capacity is 20 Mbps, enough to support the traffic load) and 
the 802.11 physical rate for UE ui is fixed at 6Mbps while that 
of UE U 2 is varied. Using ( |44| ) to account for the cost of using 
the LTE connection we expect the traffic to be offloaded to the 
802.11 network as long as its capacity is sufficient to meet the 
load. Eigure [^illustrates the LTE air time for UEs ui and U 2 
vs the 802.11 rate available to UE U 2 - Since the 802.11 rate 
for Ml is fixed at 6Mbps, mi must use the LTE connection in 
order to meet its traffc demand and indeed, as expected, it can 
be seen in Eig 11 a) that ^ and i (the fractions of LTE 


airtime used by UE mi on each subband) are always non-zero. 
In contrast, it can be seen from Eig[TTJb) that the LTE airtime 
used by UE U 2 falls to zero once the available 802.11 rate 
rises above 10Mbps i.e. once sufficient 802.11 bandwidth is 
available to support the traffic load it is fully offloaded from 
LTE onto the 802.11 link. Observe also that the 802.11 usage 
by UE Ml falls to zero when the 802.11 rate is 10 Mbps. This 
allows UE M 2 to make full use of the 802.11 network to meet 
its traffic demand, without incurring the overhead of collisions. 


B. 802.11 Multihoming 

We next illustrate how our framework can be used to 
manage 802.11 multihoming. Consider an example with one 
BS, two APs (oi and 02 ) and two UEs (mi and M 2 ). The 
available data rates are summarised in Table |Vl| Each user 
can individually achieve 10Mbps over the LTE network (and 
5Mbps each when sharing the LTE network). The users have 
a physical rate of 36Mbps and 54Mbps over the first 802.11 
WLAN and of 36Mbps and 18Mbps over the second WLAN, 
so user M 2 has a better connection via AP ai and user mi will 
have a better connection via AP oi by avoiding contention with 
user M 2 . Note that although 54Mbps link has a higher physical 
rate compared to that of the 36Mbps link, when the duration 
of a successful transmission is held constant the achievable 
throughputs are roughly equal due to the MAC framing and 
contention overhead. 


Table VII details the rate allocation found using update 
( |4 T] i. It can be seen that UE ui is mainly associated with 
AP 02 iza 2 ,ui = 0.99) and M 2 with AP oi (Zai,U 2 = 0.99), as 
might be expected in view of the link characteristics. The LTE 
network is shared equally by both users (Cb.ui = 0.5 = Ch,« 2 )' 

EigurefT^ shows how the allocation changes as the capacity 
of the 802.11 link between AP 02 and UE M 2 varies. It can 
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be seen that as the rate on the 02 -M 2 link increases UE M 2 
increasingly makes use of this link (Za 2,2 increases while Zai ,2 
falls) and conversely traffic for UE mi increasingly makes 
use of AP oi (Zai,i increases while ^ 02,1 falls) and of the 
LTE network ^ and both rise). When both UEs have 
equal rates of 36Mbps and 54Mbps through APs oi and 02 
respectively, they make use of APs equally 




(a) UE ui (b) UE U2 

Pig. 12: Illustrating 802.11 multi-homing as the rate on the 
UE M 2 - AP 02 802.11 link is varied. 


C. LTE Multihoming 

We next consider an example with two BSs (61 and 62 ), 
one AP, and 4 UEs (mi, M 2, M3, and M4). The physical rates 
available on the 802.11 and LTE links are summarised in Table 
ra We also consider the conventional user assignments 
to LTE BSs i.e. each UE associates with the LTE BS that 
provides the UE with the maximum received power. The 
achievable LTE rates when using this maximum received 
power assignment are also summarised in Table |VIII| 

TABLE VIII: LTE multihoming example: data rates. 



PHY Rales 
[Mbps] 

Technology Rales 
[Mbps] 



BS 

AP 

Lit. £>i only 

Li b 02 only 

Lib (Maximum Rx Power) 

802.11 only 

Opiimised Multi-RAT 

Ml 

26 

25 

54 

6.5 

6.25 

8.27 

7J1 

15.36 


10 

25 

27 

2.5 

6.25 

2.4 

3.75 

12.5 

113 


29 

54 

1.25 

7.25 

3.56 

7.51 

14.8 

U4 

11 

10 

13.5 

2.75 

2.5 

2.29 

1.88 

5.5 


Figure [T^ shows the proportional fair LTE sub-band al¬ 
locations and 802.11 association probabilities obtained using 
update ( [4 T] i. From Table |VIII it can be seen that a relatively 
large physical rate of 54 Mbps is available to UEs mi and M 3 
and so it can be seen from Figure [T^ that in the proportional 
fair allocation these UEs make use of the 802.11 network 
but not the LTE network, so freeing up capacity in the LTE 
network for UEs M2 and M4. UEs M2 and M4 make use of the 
LTE rather than the 802.11 network (so reducing collisions 
and increasing the 802.11 capacity available to UEs mi and 
M3), and share the LTE network evenly between them. 

It can be seen from the right-hand column of Table VIII that 
all users benefit from this use of multi-homing. In addition ob¬ 
serve that a simple aggregation of both LTE (Maximum Power 
association) and 802.11 resources results in 13% reduction in 
the proportional fair rate objective compared with that of a 
near-optimal solution. 



Fig. 13: LTE multihoming, user association probabilities per 
sub-band of the LTE BS (Q „) and per 802.11 AP (Za,u) given 
by Update (HTb. 


VII. Summary AND Conclusion 

We consider proportional fair rate allocation in a hetero¬ 
geneous network with a mix of LTE and 802.11 cells which 
supports multipath and multihomed operation (simultaneous 
connection of a user device to multiple LTE BSs and 802.11 
APs). We show that the utility fair optimisation problem is 
non-convex but that a global optimum can be found by solving 
a sequence of convex optimisations in a distributed fashion. 
The result is a principled approach to offload from LTE to 
802.11 and for exploiting LTE/802.11 path diversity to meet 
user traffic demands. 


Appendix A 
A. Optimisation problem 

With the obvious abuse of notation we add subscripts to 
the functions and to streamline the presentation since 
there is no risk of confusion. Details of the objective and 
constraints functions are as follows: 

/W = - XI log(s„ -f r„) 

u&A 
= s„ 

= Pa,u - Wa,u + X 
VeUa. 

_1 \ 

(x) = X 

ueu beB 

- X X Cb ,log (1 + ^) 

beB iex 
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eP‘^’'^Ca,n, gi^l(x) = 0 


aGAu 


,(3) 


%,nM = H n 

Oa€ V&Aa. 

•p„_i(WA{«}) 

JlWa u 

3u\^) = E TT^’ 


a^Ai 


Appendix B 

B. Proof of Theorem 

Proof: Update generates a sequence {xfc,fc = 

1, 2, • • • }. By the convexity of we have that that 5 *^*^ (x) > 
g(*)(xfc) + dg^'^\xk){x - Xfe) and so Cx^ C C (for any 
X G Cxfe we have — ( 7 ^*^(x) < h^^\x) — g^^\-x.k) — 

dg^^'>{xk){x — Xk) < 0, i = I,-- - ,to). Hence, Xk G C, 
k = 1,2, • • •. Further, /(x^+i) < /(x^) and so the sequence 
{/(xfe)} is decreasing. Since / is convex (so continuous) 
and C C B is bounded then /(x^) is bounded. Hence, by 


the monotone convergence of bounded sequences (e.g. \ 16 


Theorem 16.2]), sequence {/(x^)} converges to a finite limit. 
Let /oo denote this limit and let Coo = {x G C : /(x) = f^o} 
denote the corresponding set of limit points in C. 

From we have that for every x G Coo then x G 
argminxgCx/(x) (^Ise we could find a point y G such 
that /(y) < /(x) contradicting the fact that x corresponds to 
a limit point of monotonic sequence {/(x^)}). It follows that 
X satisfies the Fritz John conditions. 


Ko{x) := Xodfix) + ^ A«(9h«(x) - a 5 W(x)) = 0 

(45) 

Ki{x) := - g^"\-K)) = 0, i = !,■■■, I (46) 

with multipliers > 0, * = 0, • • • ,1. But these are also the 
Fritz John conditions for optimisation P, and so it follows that 
X is a stationary point of P. Since / is convex it is Lipschitz 
continuous on compact set C and so /(x^) —>■ /oo implies 
that for any 5 > 0 and k > ks, 3x G Coo s.t. \\xk — x|| < S 
provided ks is sufficiently large. By assumption, df, dh^'‘^ 
and dg^^^ are continuous, so Kq, Ki are continuous. Hence, 
Kq{xi,) —>■ 0, Ki{'x.k) —?► 0, * = 1, • • • ,l and we are done. ■ 


C. Proof of Theorem 

Proof: Since C is compact, then for every open covering 
there exists a finite subcovering. Let Uxgii>c 2 C'T be a covering 
of C consisting of balls of radius r and let UTG'i''c’p"f be a 
finite subcovering. Each set T G 4'' either covers C or has at 
least one neighbour E G 4*' such that T n S ^ 0 since C is 
connected. Now expand each set T G 4'' to be a ball of radius 
2r and define a new set T obtained by taking the union of the 
expanded set T with all of its expanded neighbours. In the new 
expanded covering rp' so obtained neighbours T and S have 
at least a ball of radius r in common, so T n E has volume at 
least v{r) where v{r) is the volume of a hypersphere of radius 
r in R". Since C is connected and the sets in form a finite 


covering, between any two points x, y G C there exists a path 
traversing a sequence of at most IT*'! neighbouring sets from 
Vk'. Given xj. G T then the probability that x^+i lies in a 
neighbour S is at least eriiy{r) and so every path is traversed 
with probability at least (e? 7 i/(r))l^ L That is, starting from 
any initial condition for every optimum x* of problem P a 
set W of size no greater than i/(2r) containing x* is visited 
with probability at least I . This holds for every 

r > 0. Selecting ^ = v{2r), then for any initial condition and 
every ^ > 0, for every optimum x* of optimisation P a ball 
B^{x*) containing x* is visited with positive probability and 
the stated result now follows. ■ 
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